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ABSTRACT 

We use numerical simulations of a (480 Mpc//i) 3 volume to show that the distribution 
of peak heights in maps of the temperature fluctuations from the kinematic and ther- 
mal Sunyaev-Zeldovich effects will be highly non-Gaussian, and very different from the 
peak height distribution of a Gaussian random field. We then show that it is a good 
approximation to assume that each peak in either SZ effect is associated with one and 
only one dark matter halo. This allows us to use our knowledge of the properties of 
haloes to estimate the peak height distributions. At fixed optical depth, the distribu- 
tion of peak heights due to the kinematic effect is Gaussian, with a width which is 
approximately proportional to optical depth; the non-Gaussianity comes from sum- 
ming over a range of optical depths. The optical depth is an increasing function of halo 
mass, and the distribution of halo speeds is Gaussian, with a dispersion which is ap- 
proximately independent of halo mass. This means that observations of the kinematic 
effect can be used to put constraints on how the abundance of massive clusters evolves, 
and on the evolution of cluster velocities. The non-Gaussianity of the thermal effect, 
on the other hand, comes primarily from the fact that, on average, the effect is larger 
in more massive haloes, and the distribution of halo masses is highly non-Gaussian. 
We also show that because haloes of the same mass may have a range of density and 
velocity dispersion profiles, the relation between halo mass and the amplitude of the 
thermal effect is not deterministic, but has some scatter. 

Key words: cosmic microwave background - large scale structure of the Universe - 
galaxies: clusters: general - gravitation - methods: n-body simulations 



1 INTRODUCTION 

The inverse Compton scattering of free electrons in the intr- 
acluster plasma with the photons of the cosmic background 
radiation produce secondary fluctuations in background ra- 
diation temperature maps. The electron motion is due both 
to thermal motions within the plasma, and to coherent bulk 
motions of the plasma (e.g. Sunyaev & Zeldovich 1980). Sim- 
ulations suggest that the distribution of these temperature 
fluctuations will be quite non-Gaussian (Seljak, Burwell & 
Pen 2001), and that the best hope of detecting this effect is 
in the tails of the distribution (e.g. da Silva et al. 2000). In 
this paper, rather than studying the distribution of temper- 
ature fluctuations in random pixels, we study the distribu- 
tion of fluctuations in pixels which are peaks. We show that 
the distribution of these peak heights should be highly non- 
Gaussian, and that the exact shape can be computed if the 
spectrum of initial density fluctuations is known. Compari- 



son with estimates of the peak height distribution from large 
n-body simulations shows that our model is quite accurate. 
This means that observations of the peaks alone should allow 
one to place constraints on what the initial conditions must 
have been. In this paper, we restrict attention to the simplest 
case in which signals come from a single redshift, chosen to 
be close to the present epoch, when nonlinear structures 
have already grown. Work in progress studies contributions 
from higher redshifts; even when a range of redshifts con- 
tribute, our conclusion that the final distribution of fluc- 
tuations should be non-Gaussian, and related to the initial 
density fluctuation distribution, should remain unchanged. 

This paper is organized as follows. Section ^ describes 
our simulations briefly. Section ^| studies the kinematic ef- 
fect, and Section ^ studies the thermal effect. The thermal 
effect is closely related to the optical depth, and so this sec- 
tion also includes a study of the optical depth distribution. 
A final section summarizes our findings. 
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Figure 1. Scatter plot showing the height of a ksz peak as a 
function of optical depth r. 



2 THE SIMULATIONS 

The simulation we will use was recently carried out by the 
Virgo Consortium. It uses 512 3 particles in a cosmologi- 
cal box of 480/t _1 Mpc on a side. The cosmological model 
is flat with matter density Qo = 0.3, cosmological con- 
stant Qa = 0.7 and expansion rate at the present time 
Ho = 70km~ 1 Mpc~ 1 . It has a CDM initial power spectrum 
computed by CMBFAST (Seljak & Zaldarriaga 1996), nor- 
malized to the present abundance of galaxy clusters so that 
<78=0.9. We emphasize that, if the kinematic effect arises 
from the peculiar velocities of clusters, then a large sim- 
ulation box such as ours is essential for studying this ef- 
fect: smaller boxes miss a significant fraction of the power 
which generates velocities, so they are liable to underesti- 
mate the magnitude of the effect (see Sheth & Diaferio 2001 
for more discussion on the finite box size effect). In addi- 
tion, our large simulation contains a large number of massive 
haloes which smaller simulation boxes are likely to misesti- 
mate. The proper population of massive halos is essential for 
studying the distribution of the thermal effect (e.g. Refregier 
& Teyssier 2000). 

From the simulation outputs we create maps of the 
thermal SZ effect, the kinematic SZ effect, and the Thom- 
son optical depth. We compute the local gas density and 
velocity from those of dark matter, and the gas tempera- 
ture is computed from the local dark matter velocity disper- 
sion, following the procedure outlined in Diaferio, Sunyaev, 
& Nusser (2000). We project the simulation box on a fine 
2400 2 grid in a random direction; thus, each pixel in the grid 
is L = 200h~ kpc on a side. This pixel scale L is about an 
order of magnitude larger than the 30/i _1 kpc spatial resolu- 
tion of the simulation. 




Figure 2. Distribution of ksz peak heights at fixed optical depth. 
Curves show Gaussians with the same rms. 



3 THE KINEMATIC EFFECT 

The optical depth is defined by 

f n\ j; ~ CTT ^ b f Jl PdmW /1 v 

T = (TT n e lectron(i) dl « TpPcrit / dl , (1) 

J flproton "0 J Pcrit 

where or is the Thomson scattering cross-section, 
w c icctron(0 is the density of electrons along the line of sight 
at I, pdm(0 is the mass density along the line of sight at I, 
Pcritfio is the average mass density of the background, and 
ilb denotes the abundance of baryons. In what follows we 
have set fib = 0.0125/ft 2 . The second equality shows the 
approximations used by Diaferio, Sunyaev & Nusser (2000) 
to relate the properties of the dark matter particles in their 
simulations to those of the gas: the gas traces the dark mat- 
ter. 

The kinematic effect temperature fluctuation ksz is de- 
fined by 

AT fib °T f Pdm(Q v(l) 

ksz = — = pcrit / dl , (2) 

1 \L(_) m pI -oton J Pcrit CligM 

where w(i)/cu g ht is the bulk velocity along the line of sight at 
I, divided by the speed of light. Notice that if one thinks of 
the optical depth as a measure of the density, then the kine- 
matic effect is a measure of the momentum. Below, we first 
describe what is seen in simulations of the kinematic effect, 
and then we present a model which provides a reasonably 
good fit to the simulations. 

Diaferio et al. (2000) presented a study of peaks in their 
simulated maps of the thermal and kinematic SZ effects. The 
panels on the left of their Figure 3 show that the mean peak 
height for the thermal SZ effect scales with the square of 
t, the optical depth: t S z = AT(K)/T(K) = -2r 2 . We will 
consider this in the next section. This section is devoted to 
a study of the relation presented in the panels on the right 
of their Figure 3. 
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Figure 3. The rms value of fcgz at fixed r as a function of r. 
Dashed line shows the scaling expected if c(fcgz|r) increases lin- 
early with optical depth. 

Fig. |l] shows our version of their Figure 3: it shows the 
heights of the peaks in the kinematic effect, fcsz, as a func- 
tion of optical depth. It was constructed following the proce- 
dures Diaferio et al. outlined; the only difference is that here 
we use a ~40 times larger simulation described in the previ- 
ous section. Notice that the mean fcsz peak height is zero for 
all t, but the scatter around zero increases as r increases. 
Fig. ti shows that, at fixed r, the distribution around the 
mean is well fit by a Gaussian. To a good approximation, 
the width of the Gaussian increases linearly with r. This is 
shown in Fig. ^. The distribution of the peak heights is got 
by integrating over all r: 

p(fcsz) = J drp(T)G(k S z\r), (3) 

where G denotes the Gaussian distribution of fcsz at fixed 
r. Because the width of the Gaussian depends strongly on 
r, the resulting p(fcsz) is a summation of Gaussians of dif- 
ferent dispersions. This is the fundamental reason why the 
distribution p(fcsz) is likely to be highly non-Gaussian. The 
histogram in Fig. W shows this distribution — note how cen- 
trally cusped it is. 

In what follows, we will describe why the distribution 
of fcsz at fixed r is Gaussian, and why o~(fcsz|i~) oc r. We 
will then provide a model for p(r). When inserted into the 
expression above, these ingredients allow us to model the 
distribution of fcsz — this produces the solid curve shown in 
Fig. g 

3.1 The model 

If we think of the integral along the line of sight I as a sum 
over discrete pixels, then we can ask how the density and 
the velocity in pixels change as we step through I. As we do 
this, it is likely that some pixels along the line of sight will 



Figure 4. Distribution of fcsz fluctuations from sources at z = 0. 
Histogram shows the measurements in the ra-body simulation, 
solid curve shows what our model predicts (eq. below). Dashed 
curves show the predicted contribution to the solid curve from 
haloes with mass in the range 10 12 — 10 13 (narrowest distribution) , 
10 13 - 10 14 , and 10 14 - 10 15 A/ Q /h (broadest distribution with 
lowest peak). 



contain, or lie within, dark matter haloes. Because the virial 
motions within haloes are random, the component of v(l) 
which is due to the virial motions will fluctuate around zero. 
The result of integrating over all I is that these internal virial 
motions will cancel out, so only the bulk peculiar velocities 
of the haloes contribute to the line of sight integration. In 
addition, there may be contributions to the integral from 
pixels which do not contain haloes. However, the density 
within a halo is on the order of two hundred times p c rit- 
So, provided that the motions of haloes are not significantly 
smaller than the motions in less dense regions, and provided 
that velocities are not correlated over scales which are on 
the order of two hundred times larger than the size of a 
typical halo, pixels associated with haloes contribute much 
more to the line of sight integral than do the much less dense 
pixels which have nothing to do with haloes. In this respect, 
our analysis is similar to that in Cole & Kaiser (1988) and 
Peebles & Juszkiewicz (1998). 

Haloes are rare, so, for sufficiently small pixels, it is 
unlikely that a given line of sight will contain more than 
one halo. Moreover, the haloes have central density cusps 
(e.g. Navarro, Frenk & White 1997). Therefore, peaks in 
the fcsz distribution occur in those pixels which contain the 
density cusps. We assume that there is one density cusp per 
halo, so the number density of peaks in the fcsz distribution 
should approximately equal the number density of haloes. 
Moreover, we can replace the integral over the line of sight 
by the value the integrand had at the position of the peak: 

fcsz=(— — J — Phalo -Rhalo = (W/Clight) T, (4) 

\Clight J "0 "Iproton 
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where v is the bulk velocity of the halo, cn g ht is the speed 
of light, phaio and 7?h a io are intended to denote the average 
density contributed by the halo over the range of pixels along 
the line of sight, -R na i , it occupied. The final expression de- 
fines t, the optical depth of the pixel containing the central 
density cusp of the halo. So, to model the distribution of 
fcsz, we must be able to model how the optical depth r and 
the velocity v depend on halo mass m. In addition, we must 
be able to estimate the number density n(m) of haloes. 

The optical depth in a square pixel of side L which con- 
tains the central density cusp of a halo of mass m and central 
concentration parameter c, can be estimated by setting 



r(m) 



dw w 
W W 



where W — [L/r s 



t(w\c, m), 



7r, and 



(5) 



t(w\ 



cm = 



0"T 



TTlproton 
0"T 



TTlproton 



-^r 3 (m)p s (m) 



dl p(x\c, m) 
2[1- 



h(w)\ 



(6) 



(e.g. Cramphorn 2001). This expression for T(w\m) is got by 
expressing all distances in units of r s (m) and then projecting 
a Navarro, Frenk & White (1997) profile along the line-of- 
sight: 



p(x\c, m) 
p s (m) 

r s (m) 

c(m) 
h(w) 



x{l + x) 2 ' 
PcTit$c{m) 
A„i 



where x 2 = I 2 + w 2 , 



Pcrit - 



c(m) 



c 3 (m) 



ln[l 
1 

c(m) 

0.13 



■ c(m)] 



3m 



c(m)/[l + c(m)]' 

1/3 



arccosh(l/w) 
Vl-w 2 ~ 
arccos(l/ui) 
Vw 2 - 1 " 



47rp C ritA n 



if w < 1 



if w > 1, 



(7) 



and A n i « 100 for the ACDM model we are considering here. 
The quantity c(m) is often called the concentration parame- 
ter of the halo. Massive haloes are less centrally concentrated 
(Navarro, Frenk & White 1997); we use the parametrization 
of this mass dependence given by Bullock et al. (2001). The 
final averaging over the circular beam of radius W can also 
be done analytically: 



r(m) 



0"T 



ln(W 2 /4)+2h(W) 



■ r s (m)p s (m) —— 

m pro ton fio W 2 

One of the dot-dashed curves in Fig. [| shows the opti- 
cal depth given by equation (^|). Notice that, for low mass 
haloes, L ^> r s , and the term in square brackets in equa- 
tion (B tends to a constant. In this limit r oc rip s (recall 



(8) 



that W 2 



which means that the optical depth is ap- 



proximately proportional to m. This is sensible, since, in this 
limit, the entire halo fits in the cell. The more massive haloes 
have larger scale radii and so they also have larger optical 
depths; at large m, r(m) is approximately proportional to 



.2/3 



(see Fig. ^|). For our purposes in this section, the im- 
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Figure 5. Optical depth r and thermal SZ effect tgz peak heights 
as functions of halo mass m. Lower dot-dashed curves show r(m) 
for profiles of the form given by Hernquist (1990) and by Navarro 
et al. (1997). Upper solid and dashed curves show -y/ |*sz I f° r 
Hernquist profiles when the circular velocity and the isotropic 
dispersion, respectively, are used to estimate the temperature. 



portant point is that the optical depth is a monotonically 
increasing function of halo mass. 

To show that this conclusion is not terribly sensitive to 
the details of the density run in the outer regions of haloes, 
the other dot-dashed curve shows the optical depth associ- 
ated with the profile shape presented in Hernquist (1990). 
The Appendix shows that, for this profile shape also, the 
relevant integrals can all be computed analytically. (We will 
describe the other two curves in the next section.) The Hern- 
quist profile scales as l/x/(l + x) 3 ; although it falls more 
steeply at large radii, it has the same small scale slope as 
the Navarro et al. profile. The fact that the optical depths 
for these two parametrizations of halo profile shapes are so 
similar shows that most of the contribution to the optical 
depth comes from the inner parts of the haloes, where the 
profiles themselves are similar. 

What about halo speeds? Sheth & Diaferio (2001) 
showed that, if the initial density fluctuation field was Gaus- 
sian, then to a good approximation, the speeds of dark mat- 
ter haloes are drawn from a Maxwellian distribution, even 
at late times. In addition, they showed that the rms halo 
speed could be computed from linear theory, and that it is 
approximately independent of halo mass. For example, at 
z — in the ACDM model we are considering, 



Vid(m) = 



430/^ 



1 + (m/2487)°- ; 



km/s, 



(9) 



where the halo mass m is expressed in units of 10 13 Mq/H. 
At z — a typical halo has m w 10 13 Mq//i, so that 
the results to follow are essentially unchanged if we set 
Vu = 387/V3 km/s and ignore the m dependence. This 
means that, to a good approximation, the distribution of 
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halo speeds along the line of sight, v, is Gaussian, and that 
this Gaussian is independent of the optical depth along the 
line of sight, r. 

Since the peaks in the kinematic effect have fcsz ~ 
t (u/ciight), where r and v are the values for the halo in 
the peak pixel, the discussion above implies that the dis- 
tribution of fcsz peak heights is really obtained by taking 
the product of two independent distributions. In particu- 
lar, because a halo's motion is independent of its mass, the 
distribution of v is independent of r, so the width of the 
distribution of ksz at fixed r should be proportional to r. 
This is consistent with our Fig. ^. Secondly, at fixed r, the 
distribution of ksz is really just the distribution of a con- 
stant times the distribution of v. Because halo motions are 
Gaussian, the distribution of fcsz at fixed r should be Gaus- 
sian; this is consistent with our Fig. | (It is worth adding 
that Sheth & Diaferio 2001 showed that the rms motions of 
haloes are higher in the dense regions. This explains why, 
in the lower panels in Figure 3 of Diaferio et al. 2000, the 
scatter in fcsz at fixed r is larger in the densest regions of 
their simulations.) 

Finally, because peaks correspond to haloes, the comov- 
ing number density of peaks equals the comoving number 
density of haloes. Let p(r\m) dr denote the probability that 
a halo of mass m has optical depth r. Then 



p(fcsz) dfcs 



dmn(m) 



^haloes 



drp(r|m) x 



p{v/ci = fcsz/r|r, m) dfcsz 



dmp(m) / drp(r\m) x 

p(v/c\ = fcsz/r) dfcsz 

T 

" p(v/ci = ksz/r) dfcsz , 
i W — ^ ■ — p(r 

dmp(m|r), 



(10) 



where we have written c\ instead of cr, g ht, and rihaioes = 
Jdrn(r). The second equality follows from writing 
n(m)/rihaioes5 = p(m), and using our assumption that the 
distribution of v is approximately independent of halo 
mass. The final expression follows from using the fact that 
p(m)p(T\m) — p(ro[r) p(t) and rearranging the order of the 
integrals. Because the integral over m leaves a function of r 
only, the final expression shows that, at fixed r, the distri- 
bution of fcsz is given by the distribution of v. Since this is 
Gaussian, our model produces a distribution of fcsz which, 
at fixed r, is Gaussian, in agreement with the simulations 
(Fig. |. 

To compute our model for the full distribution of fcsz, 
we need a model for the distribution of r at fixed m. We 
will do this in the next section. For now, note that if this 
distribution is sharply peaked about a mean value, say r(m) 
given by equation (pf), then we can replace p(r|m) with a 
delta function. In this case, 



p(fcs 



dmn(m) ^bz/^/^/o,)? /a 



(11) 



^haloes \/2n (Vld/C) )t (j7j) 

where we have explicitly shown what happens when we sub- 



stitute the Gaussian with width Vid for the distribution of 
halo velocities. 

The final integrand above is the product of the halo 
mass function and a Gaussian whose width increases as m in- 
creases. This sort of integral is just like that studied by Sheth 
& Diaferio (2001) in their model of why the peculiar velocity 
distribution of dark matter particles is not Maxwellian. In 
that case, the final distribution was got by summing Gaus- 
sians of different widths because the virial velocities did con- 
tribute to the statistic, and the rms virial velocities within 
haloes scale as m 1//3 . Here, the virial motions do not enter, 
but the Gaussians have different widths because the optical 
depth scales approximately as r(m) tx m 2//3 . 

In the next section, we will show that, in fact, p(r\m) is 
not a delta function. Nevertheless, it is reasonably sharply 
peaked, so that the delta function is not a bad approxima- 
tion. Of course, there is a small price to pay for the gain 
in simplicity. Our final expression above suggests that the 
distribution of fcsz at fixed m should be Gaussian. This is 
a direct consequence of our neglect of the scatter in r at 
fixed m. If we include the effects of the scatter (along the 
lines we describe in the next section) then the first line of 
equation (^) above shows clearly that, even at fixed m, 
the distribution of fcsz is got by summing up Gaussians of 
different dispersions, so it will be quite non-Gaussian. 

Finally, we need to specify how the number density of 
haloes depends on halo mass. To get the solid curve in Fig. ^j, 
we used the halo mass function n(m) given by Sheth & Tor- 
men (1999). The model provides a good description of the 
simulations at z — 0. The dashed curves show the predicted 
contribution to the distribution from haloes in different mass 
bins: the central spike of the distribution is from small mass 
haloes (10 12 — 10 13 Mq / h) , whereas the broad wings are en- 
tirely due to the more massive haloes (> 10 14 Mq/K). 



4 THE THERMAL EFFECT 

The CMB temperature fluctuation due to the thermal effect 
is 



tsz 



AT 
~~f~ 



fc B T c 



CTT ' 



Al 



k B T c 

2 

m< = C lieht 



dl 



me C light 
(as X — > 0) 



<7 T !T2 b f cr dm (0 . . 

-2 — / dl pdm(t), 

m c Q J 2c 2 



(12) 



where x = hv/kT, the second line shows what happens in 
the Rayleigh- Jeans limit (x — + 0, the only limit we will con- 
sider in this paper), and the third line shows the approxi- 
mation Diaferio et al. (2000) used to relate the density pdm 
and line-of-sight velocity dispersion a\ m of dark matter par- 
ticles in their simulations to the density and temperature 
of electrons. Namely, we are assuming that the gas and the 
dark matter have similar spherically symmetric density pro- 
files with isotropic velocity dispersions, so the additional as- 
sumption of hydrostatic equilibrium determines the velocity 
dispersion profile of the gas uniquely. 

Before moving on, a word on the relative amplitudes of 
the thermal and kinematic effects is in order. The first line of 
equation (^) shows that the amplitude of the thermal effect 
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Figure 6. Scatter plot showing the height of a tsz peak as a 
function of optical depth r. 



depends on frequency. The kinematic effect, on the other 
hand, is independent of frequency. In the Rayleigh- Jeans 
limit, the amplitude of the thermal effect is much larger 
than that of the kinematic effect. However, the thermal effect 
effectively vanishes at a frequency near 217GHz. Because 
the kinematic effect does not depend on frequency, the non- 
Gaussian signature due to the kinematic effect is important 
near 217GHz. 

Diaferio et al. showed that, to a good approximation, 
tsz ss — 2t 2 , where the optical depth r was defined in the 
previous section. Fig. |^ shows that this relation also de- 
scribes the thermal effect peaks in our considerably larger 
simulation box. This section describes a model for why this 
happens, which allows us to compute the distribution of tsz 
peak heights. 

4.1 The mean relation between r and tsz 

Our model for the distribution of peaks in the thermal effect 
is similar in spirit to that for the kinematic effect: each peak 
in tsz is associated with the centre of a dark matter halo, 
we ignore the effect of haloes overlapping along the line of 
sight. This means that we can replace the integral over the 
line of sight in equation ([L2|) by an integral over the density- 
weighted temperature profile of the single halo of mass m 
and concentration c which happened to be in the line of 
sight: 



tsz(c, m) 



-2r, 



or fib 



(1/ 



2w dw 
a 2 (x\c, m) 

light 



p(x\c,m), (13) 



where x 2 — I 2 +w 2 , and all distances are in units of the scale 
radius r s — r V ir/c(m), with r v i r the virial radius and c the 
concentration parameter (equation Q). As before, the inte- 



gral over w represents the beam smoothing. This shows that 
the thermal effect is proportional to the integral of the den- 
sity times the internal velocity dispersion, projected along 
the line of sight (recall that the kinematic effect is propor- 
tional to the density times the line-of-sight bulk motion of 
the halo). 

In principle, there is some freedom in deciding what 
to use for the velocity dispersion a 2 within a halo. Be- 
cause we are assuming, as Diaferio et al. (2000) did, that 
the gas and the dark matter have similar spherically sym- 
metric density profiles with isotropic orbits, the additional 
assumption of hydrostatic equilibrium determines the gas 
velocity dispersion profile, and hence the gas temperature 
profile uniquely. For the Navarro et al. (1997) profile, as well 
as for the Hernquist (1990) profile which we discussed ear- 
lier, the velocity dispersion profile a 2 (r) has a fairly compli- 
cated functional form (Hernquist 1990; Cole & Lacey 1996). 
Although the shape of <r 2 (r) is different from the circular 
velocity profile Gm(< r)/2r, the two profiles are not very 
different. The circular velocity profile has a much simpler 
functional form, and Cramphorn (2001) shows the result of 
setting <r 2 (r) = Gm{< r)/2r, and then evaluating the re- 
quired integral above for the Navarro et al. (1997) profile 
numerically. He states that the result is well approximated 
by t 2 « — tsz/2, in agreement with the simulation result in 
Diaferio et al. 

For the Hernquist (1990) profile we discussed earlier, the 
integral above can be done analytically, both for the correct 
one-dimensional dispersion, and when the circular velocity 
is used to approximate the dispersion (see Appendix A). Re- 
call that, on small scales, Hernquist's profile has the same 
slope as that of Navarro et al. (1997), so we expect that 
using it provides a reasonable analytic approximation. The 
upper solid curve in Fig. |B| shows our analytic expression for 
the beam averaged \/— tsz associated with Hernquist profiles 
when the squared circular velocity is used for a 2 , and the 
dashed curve shows the result of using the actual dispersion. 
Recall that the other curves in the Figure show the optical 
depth. Thus, Fig. [i] shows that both |tszj and r increase 
monotonically with halo mass, and that, to a good approxi- 
mation, —tsz oc t 2 . Moreover, the dashed line shows a slope 
slightly steeper than two at high masses, as it is indeed seen 
in Fig 6. The constant of proportionality is of order two, 
consistent with the simulations. The difference between the 
solid and the dashed tsz curves can be thought of approx- 
imately illustrating what might happen if our assumption 
that the gas traces the spherically symmetric dark matter 
distribution is relaxed. 



4.2 The distribution of r and tsz 

The previous subsection showed how to estimate the beam 
averaged optical depth and thermal SZ effect due to single 
haloes of mass m. If all haloes of a given mass had exactly 
the same density and velocity dispersion profiles, then we 
would be able to translate the distribution of halo masses 
into distributions of r and tsz- In fact, haloes of fixed m do 
not all have the same profile shape: the distribution of shapes 
is well parametrized by assuming that the profile is always 
of the form given by Navarro et al. (1997), but letting the 
concentration parameter c defined earlier follow a lognormal 
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Figure 7. Distribution of optical depth peaks r in the simula- 
tions (solid histograms) and in our model (smooth solid curves). 
Dashed lines show the contributions to the total from haloes with 
mass in the range 10 13 -10 14 , 10 14 -10 15 , and > lO 15 M /fc. 



distribution (Jing 2000; Bullock et al. 2001). Therefore, we 
will assume that 



p(tsz\m) dt s 



dc/c 



exp 



ln 2 [c(t S z,ra)/c(m)] 



(14) 



where the mean concentration at fixed mass, c(m), is given 
by equation (Q), and the rms scatter, a c « 0.2, is approx- 
imately independent of halo mass. The number density of 
the thermal effect peaks of height t$z is just this times the 
number of haloes of mass m, integrated over all m: 



n(tsz)dtsz= / dm n(m) p{tsz\m) disz- 



(15) 



To proceed, we need c(tsz, m). This can be got by inverting 
the tsz(c, m) relation given above. Recall that tsz(c, m) can 
be evaluated analytically if the profile has the form given by 
Hernquist (1990), so that the Jacobian dc/dtsz can also be 
given analytically. 

The number density of peaks of height r in the optical 
depth distribution can be got analogously: 



n(r) dr = 



dm n(m)p(^c{r,rn)\rn 



dc(r, m) 
d7 



dr. 



(16) 



Figs. ^ and [s] show the distributions of optical depth 
and thermal effect peak heights in the simulations. Solid 
lines show the total number density of peaks, and dashed 
lines show the contribution to the total from haloes in small 
mass ranges. In the model, for the sake of simplicity, we com- 
pute tsz (see equation (13)) using the Hernquist profile and 
its corresponding circular velocity. The dashed lines were 
computed by locating dark matter haloes using the spherical 
overdensity algorithm with density threshold 200 (details are 
in Tormen 2001, in preparation). Simulation particles which 
reside in the dark matter haloes in a given mass range were 



Figure 8. Distribution of peak heights fgz in the thermal effect 
in the simulations (solid histograms) and in our model (smooth 
solid curves). Both simulations and model assume that the gas 
traces the dark matter. Dashed lines show the contributions to 
the total from haloes with mass in the range 10 13 -10 14 (dominate 
at small t), f0 14 -10 15 , and > f0 15 M Q /h (dominate at large r). 



marked, and then the SZ effect and Thomson optical depth 
were computed by using only the marked particles. 

These distributions are bimodal; our model fits the large 
t peak for each mass range reasonably well. This suggests 
that the large r peak is due to the central cusp of a halo, 
whereas the increase at small r is primarly due to the sub- 
structure of haloes which our model does not account for. 
Note, however, that by the time we are seeing the effect of 
halo-substructure, the dominant contribution to tsz is from 
the central cusps of less massive haloes. 

The careful reader will have noticed that a scatter in 
c at fixed m leads to a scatter in optical depth at fixed 
mass r(m). This, in turn, leads to scatter in the kinematic 
effect fcsz, which is in addition to the scatter induced by 
the fact that haloes are moving with different speeds v. In 
ignoring the scatter which is due to the distribution of halo 
concentrations at fixed m, we, in effect, assumed that the 
scatter in bulk velocities v was the dominant cause of the 
scatter in fcsz at fixed m, as it is indeed. 



5 DISCUSSION 

We presented a simple model for the distribution of peak 
heights in maps of the kinematic and thermal Sunyaev- 
Zeldovich effects. In our model, which is similar in spirit 
to that proposed by Cole & Kaiser (1988), there is a one-to- 
one correspondence between a peak in the kinematic or the 
thermal effect and the presence of a massive cluster. So the 
shape of the distribution of peaks is determined by the mass 
function of clusters, and by how clusters move. If the dis- 
tribution of initial density fluctuations was Gaussian, then 
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the motions of clusters along the line of sight should be rea- 
sonably well fit by a Gaussian, so one might have expected 
the distribution of fesz peak heights to also be Gaussian. Our 
simulations showed that, in fact, the distribution of fcsz peak 
heights is highly non-Gaussian. We argued that this happens 
because the kinematic effect is proportional to the prod- 
uct of the cluster velocity and its optical depth, the optical 
depths of clusters depend strongly on cluster mass, and the 
mass range of clusters which cause peaks in the kinematic 
effect is quite large. In our model, the cluster mass function 
and the rms motions of clusters are quantities which can 
be computed if the initial spectrum of density perturbations 
is specified. Therefore, the kinematic effect can be used to 
constrain the shape of this spectrum. 

Our model also allowed us to estimate the distribution 
of peak heights in the optical depth (which is not observ- 
able) and in the thermal effect (which is). To illustrate the 
logic, we showed what one predicts for the shapes of these 
non-Gaussian distributions if gas traces the dark matter. 
Our model was able to describe the simulations quite well. 
For the optical depth and for the thermal effect, the non- 
Gaussianity was a consequence of the non- Gaussian shape 
of the halo mass function. The main aim of this paper was 
to show how our knowledge of halo speeds and profiles can 
be used to model the thermal and kinematic SZ effects. 

Real CMB maps will have the peaks of the SZ effects 
on top of the primary fluctuations. The primary fluctuations 
are Gaussian distributed with an rms, on arcminute scales, 
of about 4 x 1CP 5 . This rms is a few times larger than the 
typical thermal SZ effect, and an order of magnitude larger 
than the typical kinematic SZ effect. However, in our ACDM 
model, the SZ effect fluctuations and the primary fluctua- 
tions dominate the angular power spectrum of the CMB on 
substantially different scales, which separate at I ~ 3000 
(e.g. Springel, White & Hernquist 2001). Thus, in principle, 
an optimal high-pass filter for the arcminute scale should 
be able to isolate the peaks due to the SZ effects from the 
primary fluctuations. 

If we put our simulation box at z=0.218, then the an- 
gular size of one pixel is 1 arcminute, and the whole box 
occupies 40° x 40° square-degrees. In Fig. ^| we show a 
small 8° x 8° piece of the whole map, so as to make smaller 
angular scale structures visible. For the analysis in Fig. llOl 
however, we use the entire map. The top panel in FigureH 
shows a simulated temperature map in which the tsz effect 
fluctuations are superimposed on the primary fluctuations 
generated at the last scattering epoch. We made a two di- 
mensional random Gaussian field for the primary tempera- 
ture fluctuations. We computed the angular power spectrum 
by CMBFAST for the ACDM model we consider here. The 
bottom panel shows the effect of applying a high-pass filter 
of scale one arcminute to this map; the large scale modula- 
tions are gone, and only small scale fluctuations remain. 

Fig. |l(j shows the distribution of temperature peak 
heights in the simulated map which has both primary and 
tsz fluctuations (dashed histogram) and in the high-pass fil- 
tered map (solid histogram); i.e., the peak height distribu- 
tions in the maps shown in Fig. ^. The Gaussian-like shape 
of the dashed curve for the unfiltered map is primarily due 
to the peak height distribution of the intrinsic peaks (we 
have verified that it is similar to the analytic formula for 
peaks in two-dimensional Gaussian random fields given in 




Figure 9. Distribution of temperatures in simulated 8° x 8° 
maps of the intrinsic CMB fluctuations with the thermal SZ effect 
fluctuations superimposed. Top panel shows the full map, and 
bottom panel shows what remains after applying a high-pass filter 
of scale one arcminute. The maps shown here are part of the full 
map computed from the entire 480Mpc//i box. 

Bond & Efstathiou 1987); the tsz peaks account for most 
of the large decrement tail. The solid histogram in the left 
panel, from the filtered map, is slightly asymmetric. If we 
assume that the positive peaks are due to the intrinsic CMB 
fluctuations, whereas the decrements are either intrinsic or 
tsz peaks, then subtracting the positive side from the nega- 
tive one should leave the tsz peak distribution. The dashed 
histogram in the right panel shows what remains after do- 
ing this subtraction. It should be compared with the true 
distribution of tsz peaks which we presented earlier in the 
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Figure 10. Distribution of peak heights in simulated tempera- 
ture maps. Dashed histogram in panel on left shows the peaks 
in the unfiltered map shown (Fig. ^ top). Solid, slightly asym- 
metric histogram shows the peak heights in the high-pass filtered 
map (Fig. M bottom). Right panel compares the actual tsz peak 
decrements (solid histogram) and the peak distribution after sub- 
tracting the positive peak height distribution from the negative 
side in the filtered map (dashed histogram). Solid curve shows 
what our model predicts. The final distribution in the filtered 
map is similar to the actual tgz distribution, and is similar to 
what our model predicts. 



paper, and is now shown as the solid histogram. The smooth 
solid curve shows what our model predicts. Notice how sim- 
ilar the dashed and solid histograms are. This demonstrates 
that by suitably filtering the map, it should be possible to 
recover the true distribution of tsz peaks, and that the our 
model does a good job of describing this distribution. 

Unfortunately, this method will not work for recovering 
the distribution of the ksz peaks. To distinguish the kine- 
matic from the thermal effect, we must use their different 
spectral properties: multi-band follow-up observations of re- 
gions with deep and highly clustered tsz decrements will 
assure the presence of a supercluster where the fcsz effect is 
largest (e.g. Diaferio et al. 2000). 

In fact, the tsz peaks will also show up in clustering 
statistics. The shape of the correlation function of the peaks 
which were imprinted on the background radiation at the 
last scattering surface depends on peak height (Bond & Ef- 
stathiou 1987; Heavens & Sheth 1999; Heavens & Gupta 
2001). The presence of tsz peaks will change this depen- 
dence because, as our model shows, the two-point correlation 
function of the SZ peaks is related to the two-point correla- 
tion function of massive clusters. Since accurate analytical 
models for the clustering of clusters exist (e.g. Mo & White 
1996; Sheth & Tormen 1999), the clustering of the peaks in 
the SZ effect can be estimated rather easily, although we 
have not done so here. The clustering of clusters depends 
differently on the initial spectrum of fluctuations than the 
clustering of intrinsic peaks does; therefore, our model for 



the tsz peak distribution will be important if one wishes 
to obtain constraints on the initial fluctuation spectrum by 
studying the two-point and higher order moments of peaks 
in the microwave background on arcminute scales. 

Finally, note that our model assumes that there is a 
one-to-one correspondence between a peak in the kinematic 
or thermal effects and the presence of a massive cluster. We 
showed that, at low redshift, this is a good approximation- 
our model describes the simulations quite well. At higher 
redshifts this assumption is likely to break down for the kine- 
matic SZ effect. This is because, at higher redshifts, there are 
fewer massive haloes with sufficiently high optical depths to 
produce obvious peaks. In particular, the number density of 
massive clusters drops faster than does the typical coherent- 
flow speed, so that, at higher redshifts, an increasing fraction 
of peaks will be caused by the velocity field, rather than by 
the density field. 

Since the kinematic peaks due to coherent flows at high 
redshift have a different angular structure from those due 
to nearby clusters, an optimal filter will be able to distin- 
guish the two types of peaks. Developing a model for this 
additional effect is the subject of work in progress. 

We thank Simon White, Rashid Sunyaev, Naoshi 
Sugiyama, Atsushi Taruya and Adi Nusser for helpful com- 
ments. This research was partly supported by the NATO 
Collaborative Linkage Grant PST.CLG.976902. RKS is sup- 
ported by the DOE and NASA grant NAG 5-7092 at Fer- 
milab. 
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APPENDIX A: PROJECTIONS OF 
HERNQUIST PROFILES 

The main text uses expressions for the density and the 
density-weighted temperature profiles of dark matter haloes 
integrated along the line of sight. For the Hernquist pro- 
file, both these integrals can be done analytically (Hernquist 
f990). (Analytic results can also be obtained for many of the 
more general profiles in Zhao 1996; in the interests of brevity, 
we have not provided explicit expressions here.) 

The density at x = r/r B , where r B is a scale radius, from 
the centre of a Hernquist profile is 



p(x) = 



2tt b 3 r* ir x(l + x) 



(Al) 



where b = r s /r v i r is the scale radius in units of the virial 
radius. Reasonable agreement with Navarro et al. (f997) 
profiles of the same m and concentration c can be got by 
setting b = V2/c°- n (Sheth et al. 2001). 

The integral over the density, which is related to the 
optical depth, is 



I(z)= / p(x)dl- 



m(l + b) 2 (2 + z 2 ) h(z) - 3 



2tt b 2 r 2 ,„ 



(z 2 - l) 2 



(A2) 



where, in the integrand, x = I + z , with x, I, and z all 
in units of the scale radius 6r v i r , and h(z) is given in the 
main text. This quantity, averaged over a circular window 
of radius R, is 



Jo 



dz z 
~~ZZ 



I(z) 



m(l + b) 2 l-h(Z) 



2nb 2 rl, 



Z 2 - 1 



(A3) 



where Z = R/(br v i T ). 

There is some freedom as to how we should estimate 
the temperature. If we assume that the halo is not rotating, 
then the quantity of interest is the density- weighted line-of- 
sight velocity dispersion. If the distribution of orbits of dark 
matter is isotropic, and the gas moves like the dark matter, 
then this is 



Gm m(l + b) 4 



6 - 65z 2 + 68z 4 



24z b 



12(1 -z 2 p 



z 2 h(z)(8z (i - 28z 4 + 35z 2 - 20) 
4(1 -z 2 ) 3 

and the average over a circle is 



(A4) 



Jo 



dz£ Gm m (1 + b) A 

~ZZ [Z> ~ b^~ 2nb 2 rl r 



3 - 17Z 2 + 22Z 4 - 8Z 6 tyZ 
12(1 - Z 2 ) s 3~ 



+ h(Z) 



15Z 2 - 20Z 4 + 8Z 6 



(A5) 



12(1 - Z 2 ) 2 

If, instead, we use one half of the circular velocity squared, 



Gm(<r) Gm(l + 6) 2 



2br v 



(l + x) 2 ' 



(A6) 



2 2r 
then we need 



where 

, . , 6 + 83z 2 + 16z 4 5z 2 (4 + 3z 2 ) , , , 

^ = 12 4 ^' 

This, averaged over a circle of radius Z is 



' dz z 
~Z~Z 



V(z) 



Gm m(l + b) 4 



br viT 2-Kb 2 r\ lY 

-3 - 14Z 2 + 2Z 
\2{Z 2 - l) 3 



15Z 2 h{Z) 
\2(Z 2 - l) 3 



(A8) 
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